The structure of fluids with impurities 
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The influence of dilute impurities on the structure of a fluid solvent is investigated theoretically. 
General arguments, which do not rely on particular models, are used to derive an extension of 
the Ornstein-Zernike form for the solvent structure factor at small scattering vectors. It is shown 
that dilute impurities can influence the solvent structure only if they are composed of ions with 
significantly different sizes. Non-ionic impurities or ions of similar size are shown to not alter the 
solvent structure qualitatively. This picture is compatible with available experimental data. The 
derived form of the solvent structure factor is expected to be useful to infer information on the 
impurity-solvent interactions form measured scattering intensities. 



Dedicated to Matthias Ballauff on the occasion of his 
60 th birthday. 



I. INTRODUCTION 

Scattering techniques are powerful approaches to in- 
fer structural properties of soft materials. In var- 
ious combined experimental and theoretical studies, 
Matthias Ballauff has contributed much insight into this 
topic by successfully comparing partial structure fac- 
tors calculated from Ornstein-Zernike equations with 
experimentally determined scattering data of dissolved 
dendrimer&i^— , bottle-brush polymers^, polyelectrolyte 
brushes 7,8 , plate- like colloids^, and nanoparticle a 10 ' 11 . In 
these earlier studies the Ornstein-Zernike equations have 
been employed in the study of pair correlations of the 
dissolved particles. However, due to complex preparation 
procedures, many fluids contain impurities, which poten- 
tially influence the solvent structure factor even in the 
case of highly diluted impurities^— . In the following 
we address the problem of how impurities can alter the 
form of the solvent structure factor. By means of general, 
model independent arguments based on Ornstein-Zernike 
equations we derive a general asymptotic form of the sol- 
vent structure factor. This form of the solvent structure 
factor and the insight gained from it is expected to be 
useful also for the analysis of future studies involving im- 
purities in fluid systems. 

In the next Sec. |Il] the setting is defined and a detailed 
mathematical derivation of the asymptotic structure fac- 
tor is given. The consequences of this asymptotic form 
are discussed in Sec. Mil Finally, conclusions and a sum- 
mary are given in Sec. IIV1 



II. MATHEMATICAL PART 
A. Setting 

Consider a spatially uniform fluid mixture of three 
components of which component 1 is referred to as the 



"solvent" whereas components 2 and 3 are called "im- 
purities". The interaction potentials Uij(r) between 
two particles of components i and j are assumed to 
be isotropic and vanishing at infinite particle separation 
(r — > oo). Note that no particular asymptotic decay of 
Uij(r — > oo) is assumed. The mole fraction Xi of compo- 
nent i €E {1,2,3} may be expressed in terms of the mole 
fraction x = xi + x-$ of impurities and the composition 
4> € [0, 1] of the impurities such that x\ = 1 — x, X2 = 4>x 
and a?3 = (1 — <f>)x. For ionic impurities with valencies 
Z2 , Z3 > of the respective components the constraint 
Z2X2 — Z3X3 of local charge neutrality of the bulk fluid 
leads to the composition (f> — 2:3/(2:2 + £3), whereas the 
composition (f> of non-ionic impurities is not restricted to 
a particular value. 



B. Solvent structure factor 

In general the partial structure factors Sij(q) with 
i,j G {1, 2, 3} of a fluid with total number density g are 
given byi£ 



where 



Sij(q) = ^/XiXjSij + XiXjQhij(q), 



An f 

h ij(<l) = — jX rhij(r)sin(qr) 




(1) 



(2) 



is the (three-dimensional) Fourier transform of the to- 
tal correlation function hij(r) = 9ij(r) — 1. Here q de- 
notes the magnitude of the scattering vector and gij (r) 
is the radial distribution function. As the name im- 
plies, the total correlation function hij(r) gives a com- 
plete information about all the correlations between a 
pair of particles of components i and j. By considering 
Sij(q) and hij(q) as the (i, j)-components of the matrices 
S_(q) and h(q), respectively, and with the diagonal ma- 
trix x :— diag(xi, x<z, x§) one can rewrite Eq. ([T]) in the 
matrix form 



S(q) = y/x(l+ gVxh(q)Vx)Vx- 



(3) 



2 



The Ornstein-Zernike equation links hij(q) to the 
Fourier transform %,j(q) of the direct correlation func- 
tion dj{r)r^: 

hj(q) =Cij(q) + g^x k Cik(q)h k j(q) (4) 

k 

or in an equivalent matrix form 

(1 + g^xh(q)^/x) (1 - g^xc(q) y/x) = 1. (5) 

Hence the matrix of partial structure factors S_(g) is given 
by 



S{q) = y/x(l - Qy/xc{q)y/x) 1 Vll- 



(6) 



Introducing 7y := g^jxixjcij (q) for given q one obtains 
from Eq. (J6j) by means of Cramer's rule the solvent partial 
structure factor 



Su{q) 



X\ 



712(1 - 733) + 713(1 - 722) + 2712713723 



(7) 



1 - 7n 



(1 - 722)(1 - 733) - 723 



C. Direct correlation functions 

Since the Ornstein-Zernike equation is a definition of 
the direct correlation function c^- (r) , it does not yield a 
method of calculating this correlation function. In order 
to proceed, the asymptotic behaviour Cij(q — > 0) of the 
Fourier transform of the direct correlation function cy (r) 
is required. It is known— that dj(r > Rij) ~ —/3Uij(r) 
with 1/(3 — k B T for sufficiently large Rij > 0. Hence 
Ctj(q) = cfj(q) + c>(q) with 



35(9) 



47T 



r rcij(r) sin(qr) 



being an even entire function of q and 



4tt 



Gy(tf) = / f rpUij(r)$m(qr). 



If (3Uij(r) = Aij/r one obtains 



47rA 



(8) 



(9) 



(10) 



with an even entire function f^(q) of g. However, if 
(3Uij(r) decays faster than 0(l/r), including power laws 
cx l/r n ,n > 2, and all short-ranged interaction poten- 
tials, c^(g) is an even entire function of q. Hence the 
Fourier transform (q) is of the general form 



c l3 (q) 



4ttA; 



(ii) 



with an even entire function /y (g) of q and with a coef- 
ficient Aij, which vanishes if /3f7y(r — > 00) decays faster 
than 0(l/r). 



D. Asymptotic solvent structure factor 

For small mole fraction x <C 1 of impurities and inverse 
scattering vector magnitudes 1 j q larger than the particle 
sizes, the general form Eq. of the Fourier transform 
Cij (q) of the direct correlation function (r) leads to the 
following expressions of 7^ defined in Subsec. Ill Bl 

711 + c^y ) 



712 ~ Qy/X<p(C$ +C$q 2 



(12) 
(13) 

7i3^ev / ^W)(4 3 ) +4 1 3 ) '? 2 ) (14) 

722 -^(-^fi +C ^ + 42 ) ? 2 ) (15) 
733 -^l-0)(-^ + 4 O 3 ) +4 3 V) (16) 



723 - QX 



where and cVV are the expansion coefficients of the 
entire parts fij(q) of %(g) (see Eq. ([IT])) in q = and 
-4-22 = IbZ%, A 33 = l B z\, A 23 = -l B z 2 z 3 , where £ B de- 
notes the Bjerrum length. In the case of non-ionic im- 
purities, i.e., Z2 — 23 = 0, no pole occurs at q = in 
Eqs. C5)-CJJ. 

In the absence of impurities, i.e., for x = 0, one ob- 
tains the well-known Ornstein-Zernike structure factor 

Su(q) = (1 - q(c$ + 4iV)r 1 from E q s - © and CCD- 
From the relation Sii(0) > 0, which holds due to the 
fact that 6*11(0) is the ratio between the compressibility 
of the solvent and that of a gas of non-interacting par- 
ticles of the same temperature and density^, one infers 
1 — QCyt > 0. Moreover, 4V < is required if the bulk 
phase diagram of the pure solvent exhibits critical points 
and if the pure solvent is uniform throughout the whole 
one-phase region of the bulk phase diagram, i.e., Sn(q) 
for q > does not diverge even close to critical points. 
Note that these conditions are not fulfilled by a fluid of 
hard spheres, because there is no critical point in the cor- 
responding bulk phase diagram, which is independent of 
temperature and exhibits only a first-order phase transi- 
tion. 

Inserting Eqs. (|12l) - pT|) into Eq. ([7]) and considering 
only the dominant contributions in the mole fraction x <C 
1 of impurities leads to the asymptotic form 



4tt£ b Z2Z 3 (o) (i) 2 



& J +<&V, (17) 



Sn(q) 



Sii(0) 



i+(e«r 1 



1 + (q/nf 



(18) 



where 



1 - q(1 - x)$ - Q 2x(^ + (1 - 0)4° 3 ) )2 

(19) 

is related to the compressibility of the solvent^, £ = 
f?4^Sn(0) i s the bulk correlation length and k = 



3 



\/8tt£bI is the inverse Debye length with I = qx(z^§ ■+ 
^3 (1 — 4>))/2 denoting the ionic strength. The quantity 



f (0) _ (0) )2 
\ c 12 L 13 ) 

'4Tr£Bc[^(z2 + z 3 f 



(20) 



measures the contrast of impurity-solvent interac- 
tions. Hence the normalized solvent structure factor 
Su(q)/ 5n(0) depends on the two dimensionless param- 
eters k£ and g 2 which are related to material parameters 
such as q, x, <f>, Z2, Z3, and Ib- The form of Eq. (ITS)) has 
previously been derived in Refi-£ within a specific model 
of electrolyte solutions. The derivation presented here 
demonstrates that the asymptotic form Eq. (p~8|) of the 
solvent structure factor S\\{q) is independent of any spe- 
cific model. The consequences of Eqs. (fT5|) - ([2T)]) and the 
physical meaning of the quantity g 2 are discussed in the 
following section. 



III. DISCUSSION 



Forms of the solvent structure factor 



Considering Sn(q) in Eq. (I18[) one immediately recog- 
nizes 

Su(0) 



Su{q) 



l-3 2 «) 2 + tez) 2 
Sii(0) 



q > k 



(21) 



l + (l-g 2 )m 2 

In particular, for non-ionic impurities, i.e., for k = 0, 
the solvent structure factor Su(q) ~ 5n(0)/(l + (£q) 2 ) 
is not altered to leading order in the mole fraction x. 
However, ionic impurities due to, e.g., alkali halides of 
ionic strength / w 1 mM in water at room temperature, 
i.e., x « 4 ■ 10~ 5 , lead to a Debye length 1/k ~ 10 nm 
which is much larger than the particle size. 

For ionic impurities it is readily seen that Sn(q) is 
monotonically decreasing with q if \g\ < 1 (see Figs.QJa) 
and (b)). If \g\ > 1 a maximum of Su(q) occurs at 
q = 9max = Ky/\g\ - 1. In this latter case Sn(q max ) = 
5 n (0)/(l - «(| 5 | - I)) 2 ) is finite for £ < ( K (\g\ - l))" 1 
(see Fig.QJc)), whereas Su(q) has a pole for £ > (n(\g\ — 
1)) _1 (see Fig. Efd)). Obviously, whenever \g\ > 1, a 
divergence of Su(q) occurs upon approaching a criti- 
cal point, where the bulk correlation length £ diverges. 
Hence, by adding ionic impurities the solvent remains 
uniform within the whole one-phase region of the bulk 
phase diagram if and only if \g\ < 1, i.e., if the con- 
trast of impurity-solvent interactions is sufficiently small. 
This can be understood in such a way that a contrast of 
impurity-solvent interactions (|.g| > 0) promotes density 
fluctuations of a particular finite wave length 27r/g max 
which occur upon approaching a critical point. The 
question on the origin of a sufficiently strong contrast 
of impurity-solvent interactions is addressed in the fol- 
lowing Subsec. IIII Bl 
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FIG. 1. Influence of dilute ionic impurities on the solvent 
structure factor Sn (q). The bulk correlation length is denoted 
by £, whereas the inverse Debye length is given by k. The 
contrast of impurity-solvent interactions is measured by the 
parameter \g\ (see Eq. ([20])). For \g\ < 1 (panels (a) and 
(b)) a monotonic structure factor is found, whereas a non- 
monotonic behaviour with a maximum close to the Debye 
length 1/k occurs for \g\ > 1 (panels (c) and (d)). In the 
latter situation the maximum diverges at a sufficiently large 
but finite bulk correlation length (panel (d)). 



B. Contrast of impurity-solvent interactions 



It has been shown in the previous Subsec. IIII Al that 
deviations from a monotonic Ornstein-Zernike-like sol- 
vent structure factor Sn(q) occur only for ionic impu- 
rities and for a sufficiently large contrast of impurity- 
solvent interactions (\g\ > 1). Non-monotonic sol- 
vent structure factors have indeed been found in heavy 
water+3-methylpyridinc mixtures in the presence of 
sodium tetraphenylborate (NaBPl^) by means of small- 
angle neutron scattering^. However, the same solvent 
under the same conditions with simple salt impurities 
(LiCl, NaCl, NaBr, KC1 and MgS0 4 ) led to monotonic 
solvent structure factor o 12 ' 13 . This experimental ob- 
servation implies that the valency of impurity ions are 
not expected to be the origin for impurity-induced non- 
uniformities of the solvent. In agreement with this, it 
has recently been concluded within a theoretical model^ 
that neither the valency nor differences in solubility give 
rise to a sufficiently large parameter \g\ > 1. 

Here it is proposed that \g\ > 1 may be reached due to 
steric effects, i.e., due to the property of (ionic) impurities 
being structure-makers or structure-breakers. For spher- 



ical particles cf^ oc (i?i + R2) 3 and c^ 3 ' oc (i?i + R3) 3 
can be expected, where Ri denotes the extension of a 
particle of component i. Hence, according to Eq. (|20|) . 
\g\ > 1 may be reached for impurities of significantly dif- 
ferent particle sizes such as Na + and [BPhJ^, whereas 
\g\ < 1 for particles of similar size like simple salts. 
Non-monotonic solvent structure factors can be mea- 
sured using small-angle neutron or x-ray scattering in 



(0) 



4 



view of the aforementioned typical Debye length 1/k and 
5max/fc = \/|<?| — 1. In the case of highly diluted im- 
purities the functional shape of the measured scattering 
intensity is dominated by the solvent structure factor. 

It is important to understand that it is not the va- 
lency but the long-ranged character (oc 1/r) of electro- 
static interactions of the impurities which gives rise to 
the relevance of the parameter \g\ for the formation of 
a non-uniform solvent structure. However, whether the 
solvent becomes non-uniform upon adding ionic impuri- 
ties is mainly determined by the size-differences between 
the impurity particles. 



IV. CONCLUSIONS AND SUMMARY 

By means of general, model independent properties of 
the direct correlation functions it has been shown that di- 
lute impurities may influence the structure of a solvent if 
and only if they are composed of ions with a sufficiently 



large contrast of impurity-solvent interactions. It has 
been argued that the latter feature requires significantly 
different particle sizes. Non-ionic impurities or ions of 
similar size are not expected to give rise to a qualitative 
change in the large-scale structure of the solvent. Ionic 
impurities with a sufficiently strong impurity-solvent in- 
teraction contrast lead to a non-monotonic solvent struc- 
ture factor with a maximum at a magnitude of the scat- 
tering vector close to the inverse Debye length (Fig. [T]). 
Upon approaching a critical point, this maximum height 
diverges at a finite bulk correlation length, which corre- 
sponds to the formation of a non-uniform bulk fluid. 

For fluids with dilute impurities the asymptotic form 
of the structure factor derived here (see Eq. (|18p) is 
expected to be a useful generalization of the standard 
Ornstein-Zernike form, because additional information 
on the impurity-solvent interactions can be extracted (see 
Eq. (|20p h Application in future experimental studies on 
such systems would be highly welcome in order to test the 
derived expression further. Finally, the Ornstein-Zernike 
approach can be applied to multicomponent impurities, 
i.e., i£ {1,2, 3,.., N}. 
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